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Abstract 

We will prove an analogue of Landau's necessary conditions [Necessary density conditions 
for sampling and interpolation of certain entire functions, Acta Math. 117 (1967).] for spaces 
of functions whose Hankel transform is supported in a measurable subset 5* of the positive 
semi-axis. As a special case, necessary density conditions for the existence of Fourier-Bessel 
frames are obtained. 
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! 1 Introduction 

While Fourier Series rely on the fact that {e^^^}k^z constitutes an orthogonal basis for L'^{—it, vr), 
Nonharmonic Fourier Series allow more general sets {e**''^'}fcgz. They can be nonuniform as in 
^ ■ Riesz Basis [21], perhaps even redundant as in Fourier Frames |T3]. On their "frequency side", 

, nonharmonic Fourier series provide nonuniform and redundant sampling theorems in spaces of 

\ bandlimited functions. As a consequence of Landau's necessary conditions for sampling and 

\C) ' interpolation of such functions [inilll], we know that sampling requires {tfcjfcez to be "denser 

than Z" and that interpolation requires to be "sparser than Z". The set Z is a sequence 

of both sampling and interpolation for bandlimited functions (this is known as the Whittaker- 
Shannon-Kotel'nikov sampling theorem). 

Likewise, let Jq be the Bessel function of order a > —1/2 and jn,a its n*'* zero. Several 
classical results in Fourier analysis have been extended to Fourier-Bessel series O [H [HI [20] . 
^ I The theory of Fourier-Bessel series is based on the fact that {x^ Ja(^jn^ax)}^^Q is an orthogonal 

\ basis for L^[0, 1]. Thus, the study of more general sets {x2 Ja{tnx)}^^Q leads naturally to 

"nonharmonic Fourier-Bessel sets" . Completeness properties of such sets have been investigated 
by Boas and Pollard [2], and some stability results concerning Riesz basis have been obtained 
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in [TTj. However, the problems that arise naturally in connection with frame theory and, in 
particular, questions related to frame and sampling density, have not been investigated up to 
the present date. We will address this question in the present paper, as a special case of a more 
general result, which gives Landau-type results in the context of the Hankel transform. 

Throughout this paper S is assumed to be a measurable subset of (0, oo). 

Consider the space Ba{S) of functions in L^(0, oo) such that their Hankel transform, 

/•oo 

Ha{f){x)= / fit){xt)^/^JUxt)dt, 



Jo 

is supported in S. The special case S = [0, 1] is an important example of a reproducing ker- 
nel Hilbert space with an associated sampling theorem [S]. Moreover, this reproducing kernel 
Hilbert space is strongly reminiscent of the classical Paley- Wiener space of bandlimited func- 
tions. In particular, with a view to solving an eigenvalue problem arising in the theory of random 
matrices, Tracy and Widom [19] have constructed a set of functions which play the role of the 
prolate spheroidal functions in this situation. Such functions are examples of doubly orthogonal 
functions in the sense of Stefan Bergman [1] . This automatically implies |18j that they solve the 
concentration problem 

nr 

^k4>k{x)= / 4'k{t)'R-a{t,x)dt, 
Jo 

where TZa{t, x) is the reproducing kernel of ;Bq([0, 1]). Once we know that such functions exist, it 
becomes natural to ask if the behaviour of the corresponding eigenvalues displays the "plunging 
phenomenon" which has been observed in association with the "Nyquist rates" described in 
terms of Beurling-type density theorems (see [10], [18] and the discussion in [5l Chapter 2]). We 
will see that this is indeed the case, even in the more general case of the space Ba{S), where S 
is a measurable subset of the positive semi-axis. 

The description of our results requires some terminology. A sequence A = {tn}^=Q is a set 
of sampling for Ba{S) if there exists a constant A such that, for every / G Ba{S), 

oo 



/■oo ' — ' 

A / |/(x)pdx<5]|/(t„)p. 

■^0 n=0 



Moreover, A is a set of interpolation for Ba{S) if, given any set of numbers {flnlj^o '^i^h 
^ la^l^ < oo, there exists / G Ba{S) such that 

f(tn) = an, for every t„ G A. 

We say that a sequence is separated if the distance between any two distinct points exceeds 
some positive quantity d > 0. For such sequences we can define densities which are suitable for 
analysis of functions supported in (0, oo). 

Definition 1 Let na{r) denote the number of points of A C (0,oo) to be found in [a,a + r]. 
Then the lower and the upper densities of A are given by the limits 

D~{A) = lim inf inf °^ ^ and D'^{A) = lim sup sup ^"^^-^ . 

r— !>oo a>0 r r^oo £j>q r 

Our main results read as follows. 
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Theorem 1 Let S be a measurable subset of (0, oo) and a > —1/2. // a separated set A is of 
sampling for Ba{S) , then 

D-{A) > -m{S). (1) 



Theorem 2 Let S be a bounded measurable subset of (0, oo) and a > —1/2. // the set A is of 
interpolation for Ba{S) , then 

D+{A) < -m{S). (2) 
vr 

A major technical difficulty in the proofs of the above results arises from the translation 
invariance of Definition 1, since we cannot appeal to the translation invariance of the eigenvalue 
problem which was used by Landau in [10]. For this reason, delicated estimates of operators 
involving the reproducing kernels of the space Ba{[a,a + r]) are required. 

We will also prove that the separation condition imples the existence of a constant B such 
that, for every / G Ba{S), 

Y,\9{tn)\' < B\\g\\\ 

n 

Theorem 1 can be seen from the frame theory viewpoint. A sequence of functions {ej}j£j is 
said to be a frame in a Hilbert space H if there exist positive constants A and B such that, for 
every f ^ H, 

A\\f\\l<Y,\{f^^j)\'<B\\ffH. (3) 

Accordingly, we say that {{tnx)'^ J^^tnx)} is a Fourier- Bessel frame if there exist positive con- 
stants A and B such that, for every / G Ba[{0, 1)], 

A \f{x)\^dx<y] / itnX)-2fix)MtnX)dx <B \fix)\^dx. 

Jo -^0 -^0 

By choosing S = (0, 1) in Theorem 1 one concludes that, if {{tnx)^ Ja{tnx)} is a Fourier-Bessel 
frame, then D~{A) > In particular, the orthogonal basis {{jn,ax)^ Ja{jn,ax)} is a Fourier- 
Bessel frame, since the norm of each element is bounded away from zero and infinity [51 (2.3)]. 
It is well known (as a consequence of the Paley- Wiener theorem: see [9l Theorem 2]) that 
if / € jBq[(0, 1)] then t~"~2/(t) belongs to the Paley-Wiener space. Thus, every sufficient 
condition for Fourier frames also holds in the case of Fourier-Bessel frames. For an account of 
such conditions see , for instance, those in [13l pg. 791] and the references therein. Such an 
observation may be useful in the construction of the "Bessel analogues" of the hyperbolic lattice 
in \16\ Theorem 3.4] 

Recently, Marzo |12] applied Landau's ideas to the proof of Marcinkiewicz-Zygmund in- 
equalities in the sphere. From his work we borrow an idea to start with the estimations leading 
to inequality ([T|) and a method to deal with the case where a sequence of both sampling and 
interpolation is unknown (as in our more general situation) or do not exist (the case for higher 
dimensions in [12]). It is worth noting that analogues of Landau's necessary conditions have 
been also studied [TJIB] using techniques from time-frequency analysis [H]. 

The outline of the paper is as follows. In Section 2 we collect some results about the 
convolution structure associated with the Hankel transform. The key section is Section 3, where 
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the eigenvalue problem is formulated and the estimates of the trace and norm are obtained. 
Section 4 contains the lemmas which are required, in the proofs of the main results, to establish 
the connection between the sampling and interpolation concepts and the eigenvalue problem. 
We prove our main results in Section 5. 



2 Bessel functions and their convolution structure 

In this section we will use [15] and |16j as reference sources for some definitions and properties 
that are useful in the harmonic analysis associated with the Hankel transform. For a > — 1, the 
Bessel functions are defined by the power series, 

(-ir(f)^"+" 

Bessel functions are solutions of the second order differential equation 

(Py 1 dy / Q!^\ ^ 
dx"^ X dx \ x"^ ) ^ 

The derivative of a Bessel function can be related to a Bessel function of different order via the 
formula 

1 / d \ ™ 

Ja{x) = X jQ,_m(x), (5) 



x"^ \dx 

valid for every positive integer m. We will make extensive use of the asymptotic formulae |20] : 



Ja{x) = \ — (sinr/a; + p{x)) , (6) 

V 7^X 

with p{x) = 0{x~^), where rjx = x — (^a — |) vr, and 

J'a{x) = (cos + pi{x)) (7) 

with pi{x) = 0{x~^). Sometimes it is convenient to renormalize the Bessel functions in the 
following way: 

= r(a + l) ( - ) Ja{x). 



The functions ja are the spherical Bessel functions. They satisfy ja(0) = 1 and < 1, 

for all X € (0, oo). For the Harmonic Analysis associated with the Hankel transform one defines 
a "Hankel modulation" (m>, /) (x) = ja{^x)f{x) and associates with it a "Hankel translation" 
Ha (txI) (x) = {uixHaf) (x) = ja{)^x){Haf){x). This allows to define a "Hankel convolution" 
as follows: 

Hankel convolutions are mapped in products via the formula 



x). 
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The following property of Hankel translations, which can be found, for instance, in [15], will be 
also required: if suppf? C [0,d] and r > d then 

suppr^^ C [max{0, r — d}, r + d]. (8) 

3 The eigenvalue problem 

Let 5 be a finite union of intervals and / be the interval I = [a,a + r]. 

Let D{I) be the subspace of L^(0,oo) consisting of functions supported on / and xi the 
characteristic function of /. Let Dj and Bs denote the orthogonal projections of L'^{0, oo) onto 
D{I) and Ba{S), respectively. They are given explicitly by 

Dif = XI f and Bsf = H^DsH^f 

We want to maximize, over the functions / G Ba{S), the "energy concentration" Aj given as 

, _!j\f{t)?dt 



This is a standard problem of maximizing a quadratic form and leads to the eigenvalue problem 

\k{I,S)^k{x) = BsDi<l)k. (9) 
Writing the operators explicitly and interchanging the integrals, ([9]) becomes 

\k{I,S)(t)k{x) = j^(l)kit)ws{t,x)dt, (10) 

where, for a set X, the Reproducing Kernel wxit, x) is given by 

Wxit,x) = / Ja{ts)Ja{xs)(tx)2sds. (11) 

Jx 

Multiplying both sides of (fTO|) by {xu)^ J^[xu), integrating with respect to dx in / and changing 
the order of the integrals, gives the dual problem of concentrating on S functions whose Hankel 
Transform is supported on /: 



Afc(/, S')^fc(t) = / ■ipk{x)wi{x,t)dx. (12) 
Js 

Using this duality and a change of variables gives, for /3 > 0, the identities 

XkiI,S) = Xk{S,I) (13) 
= Xki/3I,r'S). (14) 

Now set 



I ^[xJ'^{xY - xJoc{x)J'^{x) - Ja{x)J'^{x)) if t = X. 
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and observe that 



^[a,a+r](*) x) = (a + r)lZa{{a + r)t, (a + r)x) — alZa{at, ax). (16) 

We will first study the case when 5 is a finite union of intervals. Suppose S to consist of n disjoint 
intervals (61,61 + si), (6„, 6„ + s„) and write s = si + ... + Sn- As in [TO], the cornerstone 
of the proofs consists of Norm and Trace estimates of the above operators. From the above 
considerations one has 



Trace = ^\k{I , S) = j ^Wya,a+r\{x, 



x)dx 



(a + r)TZa{{a + r)x, (a + r)x) — aTZa{ax, ax)dx (17) 
y I [a + r)Tla{{a + r)x,{a + r)x) — aTla{ax,ax)dx, (18) 



i=l 



and 



Norm = Y,Xl{I,S)= f f wl^^^^.^{t,x)dtdx 

J S J s 

~ [{a + r)TZai{a + r)t, {a + r)x) — aTla{at, ax)]^ dtdx (19) 

Js Js 

= / / [{a + r)TZa{{a + r)t,{a + r)x) - ana{at,ax)f dtdx. {2Q) 

i=i j=i Jbj 

3.1 Estimation of Trace 

In this Section we will obtain the following estimation. 

Lemma 1 For a > —1/2, Trace = ^rs + 0(1). More precisely, there exists a constant L such 
that for all positive a,r,b,s we have 



vr 



< L. (21) 



Proof. We will first estimate the function 



T(n) = / na{x,x)dx = - [xJ'^{xf - xJa{x)J'^{x) - Ja{x)J'^{x)) dx (22) 
For small x, the power series expansion of the Bessel function gives 



^ ' 2"r(a + l) ^ " 



leading, for small n, to the estimate 

T(n) = - /'"o(x2"+i)da;. 
2 Jo 
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Therefore, for a > — 1/2 and small n, the integral defining T(u) is convergent. We proceed to 
estimate T(u). Using ([7]) and (jl5p one obtains, after some simplification, 



1 



7la{x,x) = — h e{x), 



vr 



(23) 



with |e(a;)| = O (^). It is possible (and it will be required for our purposes) to improve this 
estimate even more, taking into account the cancelations resulting from the changes in sign of 
e(x). Use the second order differential equation (jH) to rewrite TZa{x,x) as 



na{x,x) = ^ (^Xj'^ix)"^ -xJa{x) ^-^J^(x) - Jq(x) + ^Ja{x) ) - Ja{x)j'a{x) 
= \ (xJo,{xf + Xj'^ixf - '^Jaixf^ 

= ^{xJ^{xf + xJ'^{xf)+o(^-^^ 

^ (^xJaix)"^ + X (^Ja~i{x) - ^Ja{x) 



2\ f 1 



2 

= \xJa{xf + \xJa~l{xf' - aJa-l{xf Ja{x) + O \\ 

The third and fifth equalities in the above calculation were obtained using ^ and the fourth 
one using ([5]). Now observe that xJa{x)'^dx = TZa{u,u) = ^ + O [^), and that the same is 
true if one replaces a by a — 1. Moreover, from ^ we obtain JJ* Jq,_i(x)^ Jo(x) = 0{1). We 
conclude that 



T(n) = -u + 0{l) 



(24) 



Finally, 

Trace 



h+Si 



y / {a + r)Tla{{a + r)x, [a + r)x) — aTla{ax,ax)dx 

i=l 



(a + r)Tla{{a + r)x, (o + r)x)dx 



E 

4 = 1 

n 

E 

i=l 

n 

V ((a + r)si - as, + 0{l)) = -rs + 0(1), 

^-^ TT 



bi+Si 



alZa{ax, ax)dx 



{a+r){bi+Si) ra{bi+Si) ^ 

TZa{x,x)dx — / TZa{x,x)dx 

{a+r)bi J abi . 



entering the estimate ()24p in the fourth identity. This is the required result. 



3.2 Estimation of Norm 

This section contains the key step, which is the estimation of 

" fb,+Si fb,+ 



Norm = / / [{a + r)Tla{{a + r)t, (a + r)x) — aTZa{at,ax)]^ dtdx. (25) 

j j Jbi Jbj 
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Proposition 1 If a > —1/2, the function Norm satisfies the estimate 

Norm > —rs — Klog{r) — L. (26) 
vr 

for some constants K and L not depending on r or a. 
3.2.1 Preparation Lemmas 

We divide the technical parts of the proof of Theorem [1] in a few preparation Lemmas. 
Lemma 2 There exists a constant L such that, for any positive b, s, a, r, we have the following 

pa~\~r poo 

/ / [{b + s)na{ib + s)t,{b + s)x) -bnaibt,bx)f dtdx > -rs - L. 

Ja Jo ^ 

Proof. The result will be derived by proving that 

ra+r poo 

I I [{b + s)TZa{ib + s)t,{b + s)x)-bna{bt,bx)fdtdx (27) 

Ja Jo 

is equal to 

/ {b + s)nai{b + s)x,{b + s)x) - ana{bx,bx)dx. (28) 

J a 

and then using Lemma [T] to bound ()28p . In order to prove that ()27p is equal to ()28p we first 
notice that, if a < 6 then 

/ abTZa{ax,at)Tla{bx,bt)dt = b / -T^q, ( -z, — t ) 7^q,(z, t)(it, 

where z = bx. Since ^TZa {^z, ^t) is the reproducing kernel of Ba{0, f ), t ^TZa {^z, |t) is a 
function in /Bq,(0, |) and thus in /Bq(0, 1), since f < 1- Using the reproducing kernel property in 
So(0, 1) one gets 

/ (a a \ „ , , , a„ /a a \ , . 

/ -nc,i^-z,-tjna{z,t)dt = -na[-^z,-zj . (29) 



We just proved that, if a < 6, 

POO 

/ abTla{ax,at)TZa{bx,bt)dt = aTZa {ax,ax) (30) 
Jo 



Then, as b < b + s, expanding (j27j) and using (|30|) on each of the 3 terms, we get the desired 
equality. ■ 



Lemma 3 Let 

(■a fa+r 



and 

fa+T foo 



pa pa+r 

P{a,r) = / nl{t,x)dtdx. 

Jo Ja 

pa+T poo 

Q{a,r) = / nl{t,x)dtdx. 

J a J a+r 
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Then, there exists constants K,L,K' and L' such that, for every a and r, 

P{a,r) < Klogr + L, (31) 

and 

Q{a,r) < K'logr + L'. (32) 
Proof. Consider the following integrals 

Hi{y,u) = / nl(t,x)dtdx 

J y~u J y+u 
f-y-u f-y+u 



H2{y,u) = / lZ\{t,x)dtdx. 

Jo Jy 

ry ry+u 
I{y,u) = / TZ'^{t,x)dtdx. 

J v—u J y 



'y-u Jy 

We will prove that Hi{y,u) and H2{y,u) are both uniformly bounded for y >u and that there 
exist constants K and L such that, for every y and u with y > u, 

I{y,u) < Klogu + L. (33) 

The bound on Q is then obtained by noticing that 

g(a,r) = Hi{a + r,r) + H2ia + r,r). 

To estimate P we have to separate in cases: if a < r then 

P{a, r) < I{r, r) < K log r + L, 

and if a > r we have 

P{a, r) = H2{a, r) + I(a, r) < K log r + L. 

We will organize the estimates in two steps: the first one contains the estimates of Hi(y, u) and 
H2{y,u) and the second one of I{y,u). 

Step 1. Using formulas ([6]) and ([7]) one can assure the existence of constants K and K' such 
that 

^"(*'^)^^0^ + ^'(^^' 

for all non-negative t,x. For y < 1 the result easily follows from the estimate of the Trace. Let 
us consider y > 1. Using (|34p in the definition of Hi{y,u) we obtain contants K and K' such 
that 

H,{y,u) < K / — - dxdt + K' / (-2—2) dxdt 

Jy~uJy+u V-'' ''/ Jy—uJy+u V-^ / 



< Ku f ( ^— ^ dx + K'u [ ( ( ^ dx 

Jy+u\x-yJ Jy+u\x-yJ \x + yj 

Ku I {-^ de + K'u r (I) ( tA^^ de 



e + 2y 



2 

de 



= {K + K')j^{^ dr. 



9 



The estimate of H2{y, u) follows the same lines. 

Step 2. A change of variables in the double integral results in 

I{y,u)=l ^ y^']Zl{yt,yx)dtdx. 

^ y ^ 

Writing the integral explicitly and inserting asymptotic formulas ^ and ([7]), one sees that 

rl rl+- 

I(y^u)= / " {£.y{t,x)+£y{t,x)fdtdx, 



with 

2 sin Tjy^t cos r]yt - sin ijytx cos ijyr, 

Cy{t,x) = ^ ^ — -r—^ ^ (35) 

vr x'^ — 

and £y{t,x) = 0{y~^). As a result (and keeping in mind that u < y), I{y,u) < 21{y,u) + Lo,for 
some Lq indepedent of y and u, where 

rl rl+- 

i{y^u) = ' {£y{t,x)fdtdx. 

^ y ^ 

Writing = — {^a — ir and using (f35]) . gives 

-r f^'^y f s'm{yt + ka)xcos{yx + ka) - sm{yx + ka)tcos{yt + ka)\'^ 

i{y^ u) = I / 75 — 2 



y 



M / X sm.{y{t — x)) + {x — t) sm.{yx + ka) cos{yt + ka) ^ ^ 



t2 



X^ 



dtdx 



= J J (f+ — ) (^y sine ^— (t — x)^ — sin(j/x + /ca) cos(rt + /ca)^ dtdx, 

y 

where we are using the usual notation sinc(a;) = sin(7rx)/7rx. From the last expression it follows 
that there exists positive constants A and B, independent of y and u, such that 

/ " (^ysinc(^-(t-x))) dtdx + B < ATT^Siy,u) + B, 

^ y ^ 

where the second inequality is obtained doing a change of variables and writing 

y y+u 

S{y,u) = / / (sinc(t — x)) dtdx. 



y-" J y_ 

n n 



Yet another change of variables shows that 

f°° 2 f 2 

S{y,u) = S{u,u) < / (sinc(t — x)) dtdx + / / (sinc(t — x)) dtdx 

Jo Jo Jr- 



2 I TT I TV 2 

(sinc(t — x)) dtdx — / (sinc(i — x)) dtdx 



JO 



u u 



/ / (sinc(t — x)) dtdx, 

71" Jo Jo 
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the last equality being true because /jgsinc(t — x)'^dt = 1. Now, Landau's inequality [TTl (8)] 
gives 

u u 

f TV f TV o U / U\ 

/ / (sinc(t - x)Y dtdx > - - C log - - 
Jo Jo Vvr/ 

It follows that 

5(y,n) <Clog(^) -S. 
Thus, for some positive constants K and L not depending on r, 

Hy, u) < K log u + L. 



3.2.2 Proof of Proposition [T]. 

Proof. Since the integrand is always non-negative, the double sum in (|25p is bounded below by 
any of the single sums. Thus, 

Norm > / / [{a + r)TZa{{a + 'r)t,{a + r)x) — aR,oi{0't,ax)]^ dtdx. (36) 

J hi Jbi 

Denote each of the double integrals in the sum ([36|) by Nj. The duality (fT3]l between / and 
S gives 

ra+r ra+r 

Ni= / [{b^ + Si)na{{bi + Si)t,{bi + Si)x)-bi^Za{bit,bix)fdtdx. (37) 

J a J a 

Set 

Ma{x,t;bi,Si) = {bi + Si)Ua{{bi + Si)t, {bi + Si)x) - biTla{bit,bix). (38) 

Using Lemma [2] with b = bi and s = Si, we get that there exists a contant Lq, not depending on 
a or r, such that 



vr 



a ra+r i-a+r coo 



> -rsi - Lo - I I Ml{x,t;b^,Si)dtdx- I I Ml{x,t;b,, s,). (39) 



J a J a J a+r 



Applying the inequality (a + 6)^ < 2(a^ + 6^) to (j38|) gives 

Ml[x, t; bi, Si) < 2{bi + s.fnliibi + Si)x, {bi + Si)t) + 2bfnl{biX, bit). 
Plugging in ([39]) . and doing a change of variable, we get 

Ni > -rsi -Lo- 2P{{bi + Si)a, {h + Si)r) - 2P{bia, bir) 
vr 

-2Q{{bi + Si)a, {bi + Si)r) - 2Q(6ia, hr), (40) 

with, 

/■a ra+r 

P{a,r)= / / nl{x,t)dtdx, (41) 

Jo ./a 
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and 

Q{a,r)= / / nl{x,t)dtdx. (42) 

J a J a+r 

Therefore, Lemma [3l provides constants Ki,K2 and Li not depending on a or r such that 

N, > -rsi -Lo-Ki log((6, + Si)r) - Ki log(6ir) - K2 log((5, + Si)r) - K2 log(6ir) - Li. (43) 
vr 

Finally, from the inequality above we find constans K and L not depending on a or r such that 

Ni > -rsi - Klog{r) - L. (44) 
vr 

Plugging this on ()36p we get (for new constants K and L), 

Norm > —rs — Klog{r) — L. (45) 
vr 



4 Sampling and interpolation and the eigenvalue problem 

In this section we will prove a series of results required to connect the sampling and interpola- 
tion problem to the eigenvalue problem of the previous section. Here we will make use of the 
convolution structure associated to the Hankel transform. 

Write A < B to signify that A < CB for some constant C > 0, independent of whatever 
arguments are involved. The next proposition is the analogue of Proposition 1 in 

Proposition 2 Let S he hounded, and let he a set of interpolation for Ba{S). Then the 
points of {tn} are separated by at least some positive distance d, and the interpolation can he 
performed in a stahle way. 



Proof. From 



one has 



f{t)= [ Haf{x)Ja{tx){tx)'2dx 
JS 

/■oo roo 

\f{t)\'< \H^f{x)\Hx= \f{x)\^dx. 
Jo Jo 



Likewise, the identity 

provides a similar estimate for |/'(t)|. The rest of the proof completely follows Landau 
Now we will prove the lemmas corresponding to Lemma 1 and Lemma 2 in 



Lemma 4 Let S he hounded and {tn} a set of sampling for Ba{S), whose points are separated 
by at least 2d > . Let I be any compact set, be the set of points whose distance to I is less 
than d, and n{I~^) he the number of points of {tn} contained in Then Xni^j+){L , S) < 7 < 1, 
where 7 depends on S and {tn} but not in I. 
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Proof. To adapt the arguments of [TT], we need the convolution structure associated with 
the Hankel transform outhned in section 2. Let hhe a function with support in [0, d] such that 
we can find constants Ki , K2 such that its Hankel transform satisfies 

< {H^h){x) < 7^2^"+^, (46) 

for every x € 5. In order to construct such a function, we use the fact (see |20^ pag. 482]) that, 
if A and B are real (not both zero) and a > — 1, then the function AJa.{z) + BzJ^{z) has ah its 
zeros on the real axis, except that it has two purely imaginary ones when A/B + a < 0. Thus, 
if Zq is a complex number outside the imaginary and the real axis, then the function 

-a-l^ . X ^ ^(zo)x-("-^) J4(x) - X~"Ux)zoJUzo) 
X ^ /-CayX, Zo) — Zq 2 2 

Zq — X 

is bounded away from zero and infinity for every x € 5. Since 

n{x,t) = H^{xmJa{-t){-tf'^){x), 

then the function h defined via its Hankel transform as 



{Hah){x) = Tla{dx,zo) 
has the desired property. Now define 

/ 7T \ ~ 1 

gix) := f hix) = ' Y^^^^ f{t)T,h{t)dt. 

Since 

H^g{x) = HM *a h){x) = x-"-5(27r)ti7„/(x)i7«/i(x) 
and / € Ba{S) then clearly also g G Ba{S). It follows that 

\\9f<Y.\9{tn)\'. 
n 

Since Kix°'~^^ < Hah{x), then 

By formula ([8]), supp Txh C [max{0, x — d}, x + d] and one can write 

(2) ' i>Tl) I «*>^^'=«''V 



\\y\\ 1 



\Jx-d J 

/ /-x+d \ j-x+d 

< / r,/i(t)V"+idt / |/(t)|2dt. 

\J x—d / J x—d 



13 



Moreover, 

rx+d 



x—d 



< 

< c, 

for some constant C > 0. We have thus shown that 

fx+d 

|2 



= x2-+l||j„(x-)/?a/l(-)f 

= / x'^''+^jl{xs)Ho,h{sfds 
Js 

= C'^ f XsJa{xs)^ ( — " ^i j 
Js \ S°+2 / 

< MS) l..,.U.sf4s 

Js 

/ xsJa{xs)'^ds, for some tq. 
Jo 



\9{^)V< \f{t)Vdt. (47) 

J x—d 



Now we impose the n[I^) orthogonahty conditions: 



for every t„ € I"*". This gives g{tn) = 0, for every t„ € /^.Finally, using the separation of 
and the definition of /+, 



\f\? <\\9\?<kY, <K' \f{t)\^dt < K' / \f{t)\ 



Since K' is independent of /, the Lemma is proved. 



Lemma 5 Let S he hounded and {tn} a set of interpolation for Ba{S), whose points are sepa- 
rated hy at least d > 0. Let L he any compact set, /" he the set of points whose distance to the 
complement of L exceeds ^, and n{L~) he the numher of points of {tn} contained in L~ . Then 
Xn(i-)-i{L, S) > 5 > 0, where 5 depends on S and {tn} hut not in L. 

Proof. We have shown in Proposition 1 that the interpolation can be done in a stable way, 
thus 

\\9f<KY,\9{tn)\'. 

n 
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Now, for each ti let (pi G B(S) be the mterpolatmg function that is 1 at ti and in the rest of 
the tn, all these functions are linearly independent. Let h be the same as in the proof of Lemma 
Hand define ipi G Ba{S) by 



"a(Pl\ 



(27r)t H^h{.) J ' 

for every x £ S (recall that x~°'~2 Hah is bounded away from zero in S). By Hankel transform, 

Given / € span{'ipn }t„e/-' ^^t 5 = / *a h. Then (7 is a linear combination of (pn with tn £ I~ , 
thus g{tn) = for tn ^ I~ ■ We get, using the estimates of the proof of Lemma H] leading to 

ft„+d f- 

Wff < K\\gf <K' Wn)? <K" / \f{t)\^dt < K" / \f{t)\^dt 

so 

Once again K" does not depend on / and we are done. ■ 

The information contained in the proof of the above Lemmas, namely inequality ()47p . allows 
to show that the separation implies the upper inequality in the definitions of sampling and the 
frame properties. 

Proposition 3 Let S he hounded, and let {tn} he a set of points separated by at least some 
positive distance d. Then, there exists a constant B such that, for every f G Ba{S), we have 

J2\9itn)\'<B\\gf. 
n 

Proof. Let g G Ba{S) and h the function constructed in Lemma HI Write 

f = Ha 



(2vr)t Hah{.) J ■ 

Clearly / G Ba{S) and by construction g = f *a h. Thus the estimate ()46p gives 

Wff = WHaff < C\\{.)-^-"^Hah){.){HM.)f = C\\Hag\? = C\\g\\\ 

with C independent of g. In the proof of Lemma 4 we have seen that ii g = f *a h then (|47p 
holds, allowing us to write 

Y \9itn)\' r' \f{t)?dt = C\\ff < C'Wgf. 

„ Jx—d 
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5 Proof of the main results 



Let {tn} be a sequence with separation constant d. We will denote by the set of points whose 
distance to S is less than d and by S~ the set of points whose distance to the complement of S 
exceeds d. We will use the notation [xj to denote the largest integer smaller or equal than x. 

Using the identities (|13|) and p4p one can see that Xk-i{r,S) is also the k^^ eigenvalue of 
the problem of concentrating on the set rS the functions whose Hankel transform is supported 
in [0, 1]. The sampling theorem associated with the Hankel transform [S] (together with the O 
(2.3)]) states that {ja,n} is a sequence of both sampling and interpolation, and is known to be 
a perturbation of the set {^}- Then, there exists T = 0{1) such that, when S is the union of 

intervals, the number of these points contained in is at most [;^rm(5)J + TA^ and their 
number in S~ is at least [;irm(5')J — TA^. Now, from Lemma 5 and Lemma 6 of Section 5, 
there are 70, 5o such that 

Ali^i„(5)j+tjv('''^) < 7o<1 (48) 

\lrn.{S)]-rN-li^^S) > 60 >0. (49) 

Theorem 3 (Sampling) Let S be a finite union of intervals. If {tn} is a set of sampling for 
Ba{S), then [0, r] must contain at least {^rm{S) — ^logr — B) points of {tn}, with A and B 
constants not depending on r. 

Proof. Let is a set of sampling for Ba{S). By Lemma H] there exists 7 independent of 
r such that 

K{i)+2 < < 7 < 1- 

Prom dMD, 

As the number of eigenvalues between 6q and 7 increase at most logarithmically with r, we have 

l-r m{S)\ - TN -1- n{I) + 2 < A' logr + B' . 
vr 

Thus, 

n{I) > — rm(5) — ^logr — B, 
vr 

for some A and B not depending on r. □ ■ 

Theorem 4 (Interpolation) Let S be a finite union of intervals. If {tn} is a set of interpola- 
tion for Ba{S), then [0,r] must not contain more than {^rm.{S) — Clogr — D) points of {tn}, 
with C and D constants not depending on r. 

Proof. Let be a set of interpolation for Ba{S). By Lemma [5] there exists 5 independent 
of r such that 

A„(/)_3 > A„(/-) > (5 > 0. 

From dMD, 

\Lrra(S)\+TNir^S) < 7o < 1. 



16 



As the number of eigenvalues between 6 and 70 increase at most logarithmically with r we have 



HI) -3) 

Thus, 



^ [^r m(5) J + TN^ < C log r + D' 



n{I) < -r m{S) +Clogr + D, 
vr 

for constants C, D not depending on r. □ ■ 

To extend the result to more general sets, we can proceed as in Landau [lOl pag. 49-50]. 
In the sampling case, the result can be extended to a general measurable set S by observing 
that it suffices to prove the result for compact sets and then cover a compact S set by a finite 
collection of intervals with disjoint interiors and measure arbitrary close to the measure of 5. In 
the interpolation case, the result is extended to bounded measurable sets by approximating in 
measure from the outside by bounded open sets. 

Finally we remark that the case of functions whose Hankel transform is supported on [a, a + r] 
cannot be reduced to a case where a sequence of both sampling and interpolation is known 
to exists (since, unlike the Fourier case, our eigenvalue problem is not translation invariant). 
Nevertheless, we can still obtain asymptotic versions of the inequalities which can be used to 
prove Theorem 1 and Theorem 2. This problem is also present in |12j . who offers a solution 
which can be adapted to our setting. From Lemma 4 we know that 

#{Xj{I, 5) > 7} < n(/+) < n{I) + o(r), r ^ 00. 

and using exactly the same argument of [12l page 582], we can obtain the lower estimate 

#{Aj(/, S") > 7} > Trace -^-^ (Trace - Norm) . 

Thus, 

n(/) > —rm{S) Alogr — B — o(r) r 00. 

TT 1 ~ 7 

The estimate required for interpolation can be obtained in a similar way. Now, Theorem 1 and 
Theorem 2 are straightforward consequences of the definitions of lower and upper density. 
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